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ABSTRACT 

A planet transits an 11**^ magnitude, GIV star in the constellation Corona 
Borealis. We designate the planet XO-lb, and the star, XO-1, also known as 
CSC 02041-01657. XO-1 lacks a trigonometric distance; we estimate it to be 
200±20 pc. Of the ten stars currently known to host extrasolar transiting plan- 
ets, the star XO-1 is the most similar to the Sun in its physical characteristics: 
its radius is 1.0±0.08 R©, its mass is 1.0±0.03 Mq, V sim < 3 km s~^, and 
its metallicity [Fe/H] is 0.015±0.04. The orbital period of the planet XO-lb is 
3. 941534±0. 000027 days, one of the longer ones known. The planetary mass is 
0.90±0.07 Mj, which is marginally larger than that of other transiting planets 
with periods between 3 and 4 days. Both the planetary radius and the inclina- 
tion are functions of the spectroscopically determined stellar radius. If the stellar 
radius is 1.0±0.08 R©, then the planetary radius is 1.30±0.11 Rj and the inclina- 
tion of the orbit is 87.7±1.2°. We have demonstrated a productive international 
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collaboration between professional and amateur astronomers that was important 
to distinguishing this planet from many other similar candidates. 

Subject headings: binaries: eclipsing - planetary systems - stars: individual (GSC 
02041-01657) ~ techniques: photometric - techniques: radial velocities 

1. Introduction 

We know of five planets that transit stars brighter than V=12: HD209458b, TrES-lb, 
HD 149026b, HD 189733b, and (now) XO-lb (Charbonneau ct al. 2000 and Henry et al. 
2000; Alonso et al. 2004; Sato ct al. 2005; Bouchy et al. 2005a, and this work, respectively 
report the discoveries).^ Observations of planetary transits of fainter stars still permit a 
rather precise measure of the relative size of a star and its planet, but the brighter stars 
provide adequate flux to enable interesting follow-up observations, including observations of 
atmospheric composition (Charbonneau et al. 2002; Vidal-Madjar et al. 2004), exosphere 
extent (Vidal-Madjar et al. 2003), and infrared emission (Charbonneau et al. 2005; Deming 
et al. 2005). The five remaining systems, discovered by OGLE (e.g. Bouchy et al. 2005b), 
are too faint (V~15) for such studies. The XO project aims to find planets transiting stars 
sufficiently bright to enable interesting follow up studies (McCuUough et al. 2005).^ 

In addition to atmospheric studies of planets, transits around bright stars also provide 
the best constraints on planetary interior models, which currently cannot explain the obser- 
vations (Konacki ct al. 2005). They also can provide an empirical determination of stellar 
limb darkening and the frequency and contrast of star spots (Silva 2003). Transits provide 
simple 1-dimensional imaging of the stellar photospheres, because an inverted transit light 
curve represents a 1-dimcnsional trace through the 2-dimensional convolution of the star's 
brightness distribution with a uniform disk equal in radius to the planet. In an analogous 
manner, eclipses in principle can provide the brightness distribution on the face of the planet 
as the star acts as a knife edge cutting across the disk of the planet during ingress or egress. 
Herein, we use "transit" to describe a planet passing in front of a star, and "eclipse" to 
describe a star passing in front of a planet or another star. 



^Alonso et al (2004) referred to both the planet and its host star as "TrES-1," but for clarity we prefer 
the trailing "b" for the planets. 

^This paper includes data taken on the Halcakala summit maintained by the University of Hawaii, the 
Lowell Observatory, the McDonald Observatory of The University of Texas at Austin, and four backyard 
observatories. 
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Comparison studies will become increasingly interesting and important as different 
methods of discovering transiting extrasolar planets produce additional systems. Some 
radial-velocity surveys (Robinson et al. 2006; da Silva et al. 2006) have a selection bias 
for high metallicity stars, in order to increase the fraction of observed stars that have hot 
Jupiters (Valenti & Fischer 2005). Photometric surveys for transits are biased toward large 
depths of the transits (i.e large planets and/or small stars) and short periods, both of which 
increase detectability. 

One of the challenges of transit surveys such as the XO project is distinguishing the stars 
with transiting planets from those without. Initial discovery of candidate planets requires 
observing thousands of stars, many hundreds of times each, with ~1% or better photometry 
per epoch. A light curve consistent with a transiting planet more often is due to a grazing 
eclipsing binary star, a triple star system, or a stellar companion that is of planetary size 
but not of planetary mass (Brown 2003). 

Whereas high-resolution spectroscopy is required to measure a planetary companion's 
mass, either moderate-resolution spectroscopy or high-fidelity photometry in many cases 
can distinguish the grazing eclipsing binary stars and the triple stars (Section 4.2) from 
transiting planets. Practical considerations may determine which is the preferred method 
(spectroscopy or photometry). A spectrum will often reveal a double-lined eclipsing binary 
immediately, and a pair of spectra separated by mere hours can reveal the large amplitude 
(km s^^) radial velocity variation of a single-lined eclipsing binary (or multiple) star system. 
Such spectra can be obtained with moderate aperture (e.g. 2-m diameter) telescopes and 
with little regard for the ephemeris of the transits. On the other hand, precise photometry 
of transits demands that the candidate be observed at specific and predictable times that 
nevertheless are inconveniently infrequent for observing in a "traditional" manner by an 
astronomer who must wait months to observe in hopes that during the week or two that 
he or she is allocated on a telescope, and the sky is clear, and the sun is down, and the 
star is at low to moderate zenith angle, the planet will transit. Faced with such prospects, 
some astronomers have deployed dedicated robotic telescopes for photometric follow up (e.g. 
Sleuth of O'Donovan, Charbonneau, & Kotredes 2004; TopHat of Bakos et al. 2006a). 

We chose a different tactic: we formed a collaboration of professional and amateur as- 
tronomers. The high-quality photometry produced by amateur astronomers of the transits of 
HD 209458 and Tres-1 (Naeye 2004) immediately after their announcements by professional 
astronomers proved that amateur astronomers can produce photometry of individual stars 
that exceeds the quality obtainable by an instrument (such as XO) that was designed to 
survey thousands of stars quickly. The XO P.I. recruited four amateur astronomers (T. V., 
R. B., B. G., and P. H.) to assist with photometric follow up of candidates identified in the 
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XO data. Important criteria in selecting these individuals were 1) the observatories should 
be distributed in longitude, permitting a given transit to have a high probability of being 
observed quickly, 2) each member should have exclusive access to a CCD-equipped telescope 
and the expertise and software to produce differentially-calibrated time-series photometry, 3) 
the telescope(s) should be capable of unattended operation, in order to reduce attrition even 
though only a small fraction of targets would be genuine planets, and 4) internet access and a 
team-oriented attitude. Henden (2006) describes other amateur-professional collaborations; 
one related to transiting planets is transitsearch.org (Seagroves et al. 2003). 

Whereas some have advocated dispersing survey instruments longitudinally around the 
globe (e.g. Alonso et al. 2004), the example of XO-lb demonstrates the effectiveness of a 
survey instrument such as XO on a single excellent site such as Haleakala, especially if it is 
coupled to a longitudinally dispersed network of observatories capable of subsequent study of 
candidates individually at specific times. We hope that as the XO Extended Team (E.T.) of 
advanced amateur astronomers grows in members, its capabilities and its value to extrasolar 
planet discovery will grow geometrically. 

2. Observations 

2.1. XO Project Photometry 

McCullough et al. (2005) describe the instrumentation, operation, analysis, and pre- 
liminary results of the XO project.^ In summary, the XO observatory monitored tens of 
thousands of bright (V< 12) stars twice every ten minutes on clear nights for more than 2 
months per season of visibility for each particular star, over the period September 2003 to 
September 2005. Prom our analysis of more than 3000 observations per star, we identified 
XO-1 (Figure 1) as one of dozens of stars with light curves suggestive of a transiting planet. 
With the XO cameras on Haleakala, we observed three transits of XO-1 in 2004 and one in 
2005, on Juhan dates 2453123, 2453127, 2453143, and 2453454. From the survey photometry 
of XO-1 (Figure 2), which has a nominal standard deviation of 0.8% or 8 mmag per obser- 
vation, we determined a preliminary light curve and ephemeris, which we used to schedule 
observations of higher quality with other telescopes, as described in the next section. 



•^The analysis of XO photometry is incomplete, so a statistical analysis would be premature and would 
tend to underestimate the efficiency for discovering similar transiting planets within the data. 
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2.2. Additional Photometry 

In June and July 2005, we observed XO-1 from three observatories in North America 
and one in Europe. Telescopes of >0.2-m diameter have adequate sensitivity for these 
purposes, and a network of them is well suited to observe candidates with known positions 
and ephemerides. A few days after XO-1 was selected from the thousands of stars XO 
monitors, the XO E.T. obtained hght curves from independent telescopes sufficient to discern 
a steep ingress and egress flanking a relatively wide fiat bottom (Figure 3). With angular 
resolution ~20 times finer than the survey cameras', the E.T. observations rediiced the 
likelihood that XO-1 was a dilute triple system, i.e. one star with a fainter, eclipsing binary 
nearby, which can mimic a transiting planet (Charbonneau et al. 2004). 

On the basis of the E.T. light curves of 2005, we acquired spectra of XO-1 as described 
in the next section. After the first few spectra showed XO-lb has substellar mass, we 
obtained additional photometry (Figure 4) with the PRISM reimager (Janes et al. 2005) in 
B. V, and R filters on the 1.8-m Perkins telescope at Lowell observatory, and with the E.T. 
telescopes. The photometry from the 1.8-m telescope was acquired in a repeated exposure 
sequence of nine 5-s R images, three 10-s V images, and one 30-s B image; the R and V 
data were averaged into 0.002-day bins. The E.T. R data were averaged into 0.006-day bins. 
Binning permitted outlier rejection and empirical estimation of the noise by the scatter of 
the individual observations, which is helpful in cases such as this in which scintillation can 
be a significant and unpredictable component to noise. 

We estimate all-sky photometric B, V, Rc, and Ic magnitudes for XO-1 and several 
nearby reference stars (Figure 1 and Table 1) cahbrated using two Landolt areas (Landolt 
1992). Using a 0.35- meter telescope on photometric nights 2006 February 25 and 2006 
March 14, we measured the fluxes of 28 and 22 Landolt stars, respectively at an air mass 
similar to that for XO-1 and established the zero points of the instrumental magnitudes 
and transformation equations for the color corrections for each filter and the CCD. The Rc 
magnitudes for XO-1 and 8 reference stars differed on the two dates by a standard deviation 
of 0.005 mag. The B, V, Rc, and Ic absolute photometric accuracies are 0.04, 0.03, 0.02, 
and 0.03 mag r.m.s., including both the formal error and an estimated systematic error. 
The Tycho magnitudes for XO-1 hsted in Table 2 transform (via Table 2 of Bessel 2000) to 
Johnson V — 11.25, i.e. 0.06 mag {2-a) fainter than our estimate. 
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2.3. Spectroscopy 

In order to measure the mass ratio of the system and to determine the characteris- 
tics of the host star, we obtained spectra of XO-1 with two-dimensional cross-dispersed 
echelle spectrographs (TuU et al. 1995; Tull 1998) at the coude focus of the 2.7-m Harlan 
J. Smith (HJS) telescope and via a fiber optic on the 11-m Hobby- Eberly Telescope (HET). 
Both telescopes are located at McDonald Observatory. The HJS spectra were obtained in 
a traditionally-scheduled manner, and the HET spectra were queued. We used an iodine 
gas cell on the HET, but not on the HJS telescope. We obtained one or two high-resolution 
(A/AA ~ 60000) spectra on each of ten nights. We extracted the two-dimensional echelle 
spectra using procedures described in Hinkle et al. (2000). 

3. Analysis 
3.1. Ephemeris 

We refined our estimate of XO-lb's orbital period P from modeling eight pairs of transits 
chosen as follows: each of four transits observed by the E. T. on 2005 June 23, 2005 July 
01, 2005 July 05, and 2005 July 12 (Figure 3) was paired with two observations of the 2006 
March 14 transit (Figure 4). We solved for the best- fitting period for each pair of transits. 
The mean orbital period, P is 3.941534 days; the l-cr dispersion is 0.000027 days. The 
heliocentric Julian date of minimum light (also mid-transit) is 

tm.l. =te + PxE, (1) 

where E is an integer, and tc =2453808. 9170±0. 0011 (HJD). Identical parameters are fisted 
in Table 3 and certainly will be refined in subsequent publications. In estimating the param- 
eters of Equation 1, we used only the photometry of the eight pairs of transits mentioned 
previously. The XO survey photometry gives similar parameters but with larger uncertain- 
ties, as expected because the survey data (by design) have poorer photometric precision and 
cadence than the E. T. data. 



3.2. Radial Velocity Measurements 

We measured XO-l's radial velocities with respect to the topocentric frame using iodine 
absorption lines superposed on the HET spectra of XO-1. We modeled the extracted spectra 
using high-resolution spectra (A/AA ~ 10^) of the Sun and the Earth's atmosphere (Wallace, 
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Hinkle, & Livingston 1998) and the HEX iodine gas cell (Cochran 2000). Using an IDL^ 
implementation of Nelder and Mead's (1965) downhill simplex minimization algorithm, 
"Amoeba," we adjusted parameters of our model to fit the observations. The model includes 
convolution of our input spectra with a best-fitting Voigt profile to approximate the (slightly 
non-Gaussian) line-spread-function of the instrument. The free parameters of our model are 
a continuum normalization factor, the radial velocity of the star, the radial velocity of the 
iodine lines (which represent instrumental deviations from their expected zero velocity with 
respect to the observatory) , and an exponent (optical depth scale factor) that "morphs" the 
depths of the fines, as an arbitrary method of adjusting the solar spectrum to that of XO-1. 
Due to the iodine absorption, we could not estimate the continuum level by interpolating 
between local maxima in the spectrum, so instead we solved for the continuum iteratively, as 
required to improve the fit between our model and the observations. In the manner described 
above, for each ~15 A section of each individual spectrum within the region of the recorded 
spectrum with significant iodine absorption, 5100 — 5700 A, we estimated the radial velocity 
of the star. From the approximately normal distribution of the resulting radial velocity 
estimates for each epoch, we calculated the mean radial velocity and its uncertainty. We 
transformed our measured radial velocities to the barycentric frame of the solar system and 
phased them to the ephemeris known from the transits (Figure 5 and Table 4). 

The l-(7 internal errors of the radial velocities from the iodine spectra are ~15 m s~^ per 
epoch (i.e. per average of two 15- minute exposures) with the 11-m HET and ~65 m s~^ per 
15- or 20-minute exposure with the 2.7-m HJS telescope. We used only the iodine-calibrated 
radial velocities to determine the amplitude of the radial velocity variation, assumed to be 
sinusoidal and phased with the transits. An eccentricity approximately equal to zero is 
expected theoretically (Bodenheimer et al. 2001) and is consistent with the radial velocities 
(Figure 5c). The maximum likelihood radial velocity semi-amplitude K = 116±9 m s~^. 



3.3. Spectroscopically-Derived Stellar Properties and Planetary Mass 

An initial estimate of the spectral type of XO-1 was obtained using the fine depth ratio 
method developed by Gray & Johanson (1991). The spectral type is estimated by measuring 
the depth of the V I 6251.83 A line divided by the depth of the adjacent Fe I 6252.57 A line. 
The 4 spectra from the HJS telescope were summed together and the individual line depths 
were measured by fitting a Gaussian function to each line. The resulting line ratio is 0.14, 
which corresponds to a Gl star (Table 1 of Gray & Johanson 1991). Differences in metallicity 



^IDL is a registered trademark of Research Systems, Inc. 



-8- 



can affect tfiese results (Gray 1994), so we estimate an uncertainty of one spectral subclass 
for the spectral type of XO-1. These published line depth ratios are for dwarf stars, and Gray 
& Brown (2001) show that there is a gravity dependence to the spectral type classification 
that results from this ratio. We estimate the luminosity class V by examining the gravity- 
sensitive Na D doublet and the Mg b triplet. Below we estimate the \ogg = 4.53 ± 0.065 
cm s~^, whereas early G-type luminosity class 111 giants have logg ~ 2.5 cm s^^, which 
produces substantially narrower Na D profiles (Gray 1992). Our final spectral type estimate 
for XO-1 is GIV. 

We used the software package SME (Valcnti & Piskunov 1996) to fit each of the four 
spectra of XO-1 from the HJS telescope with synthetic spectra. We used the methodology 
of Valenti & Fischer (2005), including their minor corrections to match the Sun and remove 
abundance trends with temperature (negligible in this case) . Because of gaps between echelle 
orders, the HJS spectra are missing the wavelength intervals 6000-6123 A and 6113-6118 A, 
which were included in the Valenti & Fischer analysis. These wavelength intervals are also 
missing from our extracted HET spectra because the relevant echelle orders span the gap 
between the two detectors. 

We averaged our SME results for the four HJS spectra, obtaining the parameter values 
in Table 5. Each value in the last column of the table, labeled "Precision" because systematic 
uncertainties are not included, is the standard deviation of the four measurements divided by 
\/3 to yield the formal uncertainty in the mean. The median value of each derived parameter 
(not given) differs from the mean by less than the uncertainty in the mean. The final row 
in the table gives [Si/Fe], which Valenti & Fischer (2005) used as a proxy for alpha-element 
enrichment, when interpolating isochrones. 

Figure 6 shows XO-l's spectrum in the region of the Mg b triplet, which is the dominant 
spectroscopic constraint on gravity. These three Mg lines also have a significant impact on 
the global [M/H] parameter, which is used to scale solar abundances for all elements other 
than Na, Si, Ti, Fe, and Ni. The unidentified absorption fines at 5171 A and in the wing of 
the 6020 A line are common in stars of this type (Valenti & Fisher 2005) so presumably are 
unrelated to the presence of a hot Jupiter. 

Following the methodology of Fischer & Valenti (2005), we interpolated Yonsei-Yale 
(Y^) isochrones (Demarque et al. 2004) to determine probability distribution functions for 
the mass, radius, gravity, and age of XO-1. The trigonometric parallax of XO-1 is unknown, 
so we initially assumed distances of 180, 200, and 220 pc with an adopted uncertainty of 

10 pc in each case (which affects the width of the resulting distribution functions). We 
measured a V magnitude of 11.19 ± 0.03 for XO-1 (Table 2) so XO-1 would be at 190 pc, if 
it has the same intrinsic brightness as the Sun. 
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We used our spectroscopic effective temperature, spectroscopic gravity, and an assumed 
distance to derive a bolometric correction by interpolating tlie "liigli temperature" table 
from VandenBerg & Clem (2003). We combined the bolometric correction with the observed 
V-band magnitude to determine stellar luminosity. Then we used the stellar luminosity and 
our spectroscopic effective temperature, iron abundance, and alpha-element enrichment to 
interpolate the isochrones to produce the probability distribution functions in Figure 7. 

As expected, the most probable mass, radius, gravity, and age are essentially solar for 
an assumed distance of 200 pc. If we assume that XO-1 is at 180 pc instead, then the most 
probable age drops to 2 Gyr, which is unlikely given how slowly XO-1 rotates (V sim < 3 
km s~^). If we assume XO-1 is at a distance of 220 pc, then the gravity from stellar structure 
models deviates significantly from the gravity we derived spectroscopically. Thus, we assume 
in subsequent analysis that XO-1 is at a distance of 200 pc, but clearly the trigonometric 
parallax needs to be measured. 

Assuming a distance of 200 ± 10 pc yields a 68% confidence hmit of 1.00 ± 0.025 for 
M/M0, 1.00 ± 0.05 for R/R©, and 2.5-6.6 Gyr for the age (Figure 7). The stellar mass 
combined with our measured semi-amplitude K = 116±9 m s~^ yields a mass of 0.90±0.07 
Mj for the planet XO-lb. Within the broader range of distances of 200 ±20 pc (Table 6), the 
light curves and the spectroscopy produce consistent estimates of the stellar radius (Figure 
8). Based upon XO-l's proper motion, radial velocity, and distance, its velocity vector with 
respect to the local standard of rest is (u,v,w) = (8, 9, 23) km s^^, which is consistent with 
the kinematics of a G-dwarf with an age of a few Gyr (Wyse & Gilmore 1995). 

3.4. Light Curve Modeling and Planetary Radius 

We modeled the observed transit light curves using the algorithms of Mandel & Agol 
(2002) and adopt their notation. For the star, we interpolate quadratic limb darkening 
coefficients from Claret (2000) for the filter bandpasses: B (71 = 0.639; 72 = 0.179), V 
(71 = 0.437; 72 = 0.295), or R (71 = 0.339; 72 = 0.324), based upon our spectroscopic 
gravity and effective temperature and assuming solar metallicity, which is consistent with 
our spectroscopic determination (Table 5). In addition to the two limb darkening coefficients, 
parameters of a model light curve include the mass of the star, the size ratio p — Vp/r^ 
of the radius of the planet to the radius r* of the star, the inclination i and the angular 
frequency uj of the orbit, and the heliocentric epoch of the transit's center, tc- Note that the 
shape depends on the size ratio p — Vp/ r* and not the two radii separately. 

Of the seven input parameters (71, 72, M*, p, i, u!,tc) the most precisely measured are cu 
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and tc- As previously noted, we estimated 71 and 72 from the spectroscopically-determined 
logg and Tg//. Our derivation of the other parameters is similar to that of Brown et al. 
(2001). We begin by adopting the stellar mass and its uncertainty from the previous 
section. Combining the mass M^, and the angular frequency u with an assumption of a 
circular orbit gives the velocity of the planet passing in front of the star. Combining that 
velocity with the observed duration of the transit (from 2nd contact to 3rd contact) gives a 
chord length for the transit across the star. The half chord length L must be less than or 
equal to the radius of the star, i.e. L < r^. The stellar radius r* also has an upper bound 
from the light curve: for a transit of a given L, the larger r* becomes, the longer will be the 
ingress and the egress, as the transit becomes more grazing. Thus, the duration from 1st 
to 2nd contacts (or 3rd to 4th) creates an upper bound, r^, < l.SR©, whereas the duration 
from 2nd to 3rd contacts creates a lower bound, r* > O.QR©. 

For rhetorical purposes we have discussed the transits as if the contacts are distinct 
features in the light curves, although with the precision of our photometry, they are not. 
Instead, the model light curves become too narrow compared to the observations if < 
O.9R0, and the models become too V-shaped to match the "steep-edged" U-shape of the 
observations if > 1.3R0. In between these extremes, r^, and i are anti-correlated: increases 
in r* can be accommodated by decreases in i, and vice versa: 



where 6q = IR0/IA.U. = 4.65 x 10^^ rad, P is the planet's orbital period in years, p is the 
dimensionless mean density of the star in solar units i.e. p = (m*/M0)(r*/R0)~^, and we 
determine the scalar 5 so as to match the relation between each best-fitting inclination i for 
each assumed value for p. Geometrically, q is the dimensionless duty cycle of the transit, 
i.e. the duration of the transit divided by the period, where the transit's duration is defined 
as the interval between the two moments at which the planet's center coincides with the 
limb of the star. For XO-1, q =0.0275±0.0005, r* ^l.OiO.OSR©, =1.0±0.03 M0, so 
p = 1.0 ± 0.24 and thus i =87.7°±1.2°. 

The radius of the planet Vp is approximately proportional to r*, because the depth of the 
transit is approximately equal to p^. At a more precise level of approximation, rp oc r^, where 

K accounts for limb darkening, which necessitates an increasingly larger rp as increases 
(Equation 3). For various presumed values for r^, and m^,, we minimized of the R-band 
data by varying tc, Vp, and i. The region of high likelihood for XO-1 is described by the 
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functional relation, 

rpRi(1.3O±O.O3)(r,/r0)'^Rj, (3) 

where k ~ 1 + 0.57iL/r*. As previously defined, 71 is the linear limb darkening coefficient, 
and L is the half chord length of the transit's path across the star. For XO-1 observed in 
R band, k = 1.17. Substituting the best-fit value and uncertainty for the radius of the star 
into Equation 3, we conclude that the planet's radius^ as determined from the R-band light 
curves is 1.30±0.11 Rj. 

Limb darkening is stronger in the blue than the red, so a transit with inclination i — 90° 
should be deeper in the blue than the red. The difference in transit depths in B band 
compared to R band would be 2 mmag for XO-lb iii = 90°. At the other extreme, a grazing 
transit, the transit should be shallower in the blue than the red. The transit depth color 
difference's constraint on inclination requires more accurate photometry than that presented 
here (cf. Jha et al. 2000). 



4. Discussion 

4.1. Transit Timing 

Other planets in the XO-1 system could be revealed by their transits or by their pertur- 
bations of the radial velocity of XO-1 or on the time of arrival of the transits by XO-lb. We 
have not observed sufficient numbers of transits by XO-lb with adequate quality to justify 
an analysis of time-of-arrival perturbations. To motivate such observations, we note their 
scientific value and present an analysis of the dependence of timing precision on photometric 
precision. 

Transit timing observations would be especially sensitive to planets in orbital reso- 
nance(s) with XO-lb (Steffen & Agol 2005; Holman & Murray 2005). In principle the error 
in a measurement of the time of arrival of a transit is proportional to the photometric er- 
ror per unit time. We estimate that for an idealized transit light curve, with an isosceles 
trapezoid shape, i.e. with negligible limb darkening, the error in the time of arrival (TOA) 



^Authors unwittingly create ambiguity if they do not define Rj explicitly, because various averages of 
Jupiter's polar radius Rp and its equatorial radius Re could be selected for Rj. Because Rp ~ 0.94i?e 
(Mallama et al. 2000), such ambiguity can be significant compared to observational uncertainties; e.g. 

Knutson et al. (2006) report the radius of HD 209458 with a fractional uncertainty of 1.9%. From reported 
masses and mean densities, we determined that many authors tacitly define Rj = i?e = 71492 km; hence we 
have also, for the sake of consistency. 
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of each transit is 



a\TOA) ^ ( 



a, 



m 



W 




(4) 



where am/d is the dimensionless ratio of the standard deviation of the photometry to the 
depth of the transit, is number of observations obtained during ingress and/or egress, w 
is the time interval between first and second contacts, and e is a characteristic error due 
to (uncahbrated) secular variations in gain or atmospheric transparency over a time scale 
equal to w. Implicit in Equation 4 is that 1) the light curve shape is known, 2) there are 
sufficient observations to set the level of the light curve's out-of-transit portion precisely, 3) 
the noise has a normal distribution, and 4) N ^ 1. For example, if e = and one can obtain 
am = 2 mmag precision per 1-minute observation of a transit with depth d = 20 mmag 
and duration of ingress w = 20 minutes, then for an observation of ingress, = 20 and 
the time can be measured with a standard deviation a {TO A) = 27 seconds. If ingress and 
egress are both observed, then N — AO and a{TOA) — 19 seconds for the transit as a whole. 
For € — ±1 mmag the measured time of arrival of the transit is shifted by ±30 seconds. If 
there are additional uncahbrated systematics in an observed light curve, then Equation 4 
will underestimate the uncertainty of the time of arrival. In particular, if only ingress or 
egress is observed, then secular trends are more difficult to calibrate and the residual error 
in calibration, e, will be larger than if many observations are obtained immediately before 
and after each transit. 

Because scintillation is independent of the brightness of the star, the photometric pre- 
cision, and hence the timing precision of transits observed from the Earth's surface, is not 
strongly dependent of the brightness of the target star, for stars as bright or brighter than 
XO-1. If the star is sufficiently faint that the Poisson contribution to noise is much larger 
than the scintillation, e.g. if the sensor or the spectral filters limit the collecting efficiency, 
such as in B band, or if spectroscopy is the goal, bright stars and/or large telescopes can 
have a significant advantage. The Poisson noise for XO-1 observed with a 0.3-m telescope 
with an R-band filter is 0.15% (i.e. ~1.5 mmag) in a 60-second exposure. Near sea level, 
typical scintillation for such a telescope and exposure time is 0.9, 1.7, and 2.9 mmag at air 
masses of 1, 1.5, and 2, respectively (Young 1967; Dravins et al. 1998). For differential 
photometry, comparison stars comparable in brightness to the target star generally will be 
nearer to a faint target than a bright target. The latter effect can be important if the nearest 
suitable comparison star is too distant in angle to observe with the target simultaneously on 
the same sensor. 

We searched for a second planet transiting XO-1 by removing the observations within 
0.06 days of minimum light of the known transit. The sensitivity to transits decreases 
with increasing period, and the a priori likelihood that a planet orbiting XO-1 will transit 
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also decreases with period. Due to the measurement errors of the photometry from the 
XO cameras and the radial velocities, respectively, we would not have detected a second 
planet orbiting XO-1 if its radius Vp < 0.5 Rj or if its mass rUp, inclination i, and angular 
frequency, are such that nip sin{i) {uj2/ujif'^'^ < 0.5 Mj, where 002! ^\ is the ratio of the 
angular frequency of the second planet to that of XO-lb. The complementary statement, 
that we would have detected a second planet if either or both of those inequalities were not 
true, cannot be definitive in general. In particular, additional planets could orbit XO-1 with 
periods much longer than that of XO-lb. Monitoring XO-l's radial velocity could reveal 
additional planets and enable predictions of the approximate times that they might transit. 

4.2. A Triple Star? 

Mandushev et al. (2005) demonstrate that a triple star system (gravitationally bound or 
an asterism) can exhibit a light curve and a radial velocity curve similar to that of a planet or 
brown dwarf transiting a single star. The rotational broadening of Mandushev et al's stars, 
(V sini = 34 km s~^ for the FO primary; V sini = 15 km s~^ for the GO secondary) produced 
an oscillating asymmetry in the spectral lines that could be misinterpreted as periodicity 
in the radial velocity of a single star induced by a dark companion of substellar mass. We 
performed a simpler form of Mandushev et al's analysis for XO-1 which we summarize here. 
Our analysis is simpler because XO-1 has narrow spectral lines (V sini < 3 km s^^), so 
perturbations by a stellar companion would be observed as shifts of narrow lines, rather 
than distortions of broad lines. 

Mandushev et al's light curve consists of three parts equal in duration: ingress, flat 
bottom, egress. Although a transiting planet can produce such a light curve, the a priori 
probability of the duration of ingress (or egress) being as large or larger than the duration 
of the transit's middle section, between contacts 2 and 3, is approximately 1 — -^/l — 2rp/r*, 
which is always greater than or equal to p = r^/r*, and equals ~10% for a typical Jovian 
planet transiting a solar-type star. In practice, we expect such grazing transits to be detected 
with lower than a priori chance, because their reduced depth and duration both reduce 
detectability. On the other hand, if p > 0.375, the probability of ingress (or egress) being 
at least as long as the middle section of the light curve is > 0.5; that is, eclipsing binary 
stars typically exhibit such light curves, whereas transiting planets typically do not. Such 
probabilistic arguments cannot prove a given system is not a transiting planet, but we have 
used them to focus our efforts on "steep-edged" U-shaped light curves instead of those with 
ingress and egress comparable in length to the fiat bottom. 

For XO-l's light curve, the total duration of the fiat bottom is ~5 times the duration of 
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ingress or egress, which imphes that the ratio of diameters of the echpsing bodies must be 
> 5. For a Gl V primary (r^, = 1.0 R©), an echpsing stellar companion must have a radius 
< 0.2 Rq to match the duration of ingress relative to the flat bottom, which implies a 
spectral type no earlier than M5 V. Diluting the transit depth with a second Gl V star in 
the system would mimic the observed depth (2%) without affecting the durations.^ However, 
the radial velocity shifts of a Gl V star orbited by a O.IM© companion would be ~10 
km s~^. For a star with narrow lines hke XO-1, velocity shifts of this amphtude would be 
apparent in our spectra, in contrast to the situation of stars with much broader spectral 
hnes. 

Modeling XO-1 as a triple star with a second star fainter than the GIV primary encoun- 
ters difficulty with the radius of the third star (that eclipses the second star) . As the second 
star contributes a lesser fraction of the total hght of the system, we are forced to make the 
third star an ever larger fraction of the size of the star it eclipses, in order to match the 2% 
depth. However, from the light curve we know that the ratio of diameters of the eclipsing 
bodies must be > 5. 

Modeling XO-1 as a triple star cannot include a secondary star brighter than the ob- 
served narrow-lined GIV star, unless the secondary star's spectrum is featureless, perhaps 
due to great rotational broadening. However, the broadband colors of XO-1 are consistent 
with a slightly reddened^ Te// = 5750 K star (Table 2). 

We conclude the XO-1 transit cannot be caused by a stellar binary diluted by a third 
star that is either physically bound or along the line of sight, due to the mutual constraints 
of eclipse depth, eclipse shape, and a single-hned GIV spectrum with V sin i < 3 km s"^. 

XO-lb must be a planet, but the XO-1 system could include a wide-separation stellar 
companion that is either much fainter than or very similar in spectra to the GIV star. For 
example, HD 189733 has a faint stellar companion (Bakos et al. 2006b). For XO-1, the 
possibility of a stellar companion could be tested with infrared adaptive optics on the largest 
available telescopes. 



^In this hypothetical example, we could replace the M5V star with a planet twice the nominal mass of 
XO-lb and twice the latter's projected surface area, and all the observational constraints would be equally 
well satisfied. However, the density of the planet would be l/v^ times the nominal (single-star hypothesis) 
value. Such a contrived example seems less plausible than the nominal system of a single star and planet. 

^In the direction of XO-1, the interstellar extinction Ay = 0.04 mag as inferred from infrared observations 
with a 6' (FWHM) beam (Schlcgel, Finkbeiner, & Davis 1998). We neglect its small (<2%) effect on the 
distance and the radius of XO-1. 
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4.3. Comparison to other Transiting Planets 

Figure 9 compares the radius and mass of XO-lb with those of the nine other transiting 
extrasolar planets as tabulated by Charbonneau et al. (2006). The mean density of XO-lb 
is 0.51±0.13 g cm^^. The radius of the star XO-1 is not known to sufficient accuracy to 
assert that its planet XO-lb is truly anomalously large. If the Gl V star XO-1 has the 
same radius as the Sun, then the planet XO-lb is approximately the size of HD 209458b. 
The radius of each of these two large planets, ~1.3 times the radius of Jupiter, is ~20% 
larger than theoretical models predict (Table 1 of Bodenheimer et al. 2003). To keep the 
planet's size large, an additional heat (or power) source is hypothesized deep inside the 
planet. Knutson et al. (2006 and references therein) summarize three proposed models for 
that heat source: 1) strong winds on the planet, with the original power source being the 
incident light of the star, 2) tidal dissipation, presumably due to a second planet forcing the 
hot Jupiter into an eccentric orbit, and 3) resonances between the spin precession and the 
orbital precession, a so-called Cassini state. At 0.90 Mj, XO-lb is ~2 times more massive 
than the best- fit correlation of mass and period noted by Mazeh, Zucker, & Pont (2005), 
whereas HD 149026b is ~2 times less massive than the same correlation. 



5. Summciry 

The planet XO-lb is the first planet discovered by the XO project. It is similar to other 
transiting extrasolar planets in its physical characteristics (Section 4.3; Table 3; Figure 9). 
Of those nine others, HD 209458b has the most similar radius to the nominal radius of XO- 
lb. The host star XO-1 is very similar to the Sun in its physical characteristics (Tables 2, 

5, and 6). Equation 2 expresses the inclination of a transiting planet's orbit as a function 
of its period, the duty cycle of the transit, and the mean density of the star. For transits 
with imprecisely measured inclinations near 90°, the derived ratio of the radii of the planet 
and the star is approximately proportional to the stellar radius raised to a power equal to 
one half the linear limb darkening coefficient (Equation 3). We note the scientific potential 
of transit timing and estimate limits to its accuracy. We discuss evidence that disproves the 
hypothesis that the periodic modulations of the photometry and the radial velocity of XO-1 
are due to a triple star system. 

The University of Hawaii staff have made the operation on Maui possible; we thank espe- 
cially Bill Giebink, Les Hieda, Jake Kamibayashi, Jeff Kuhn, Haosheng Lin, Mike Maberry, 
Daniel O'Gara, Joey Perreira, Kaila Rhoden, and the director of the IFA, Rolf-Peter Ku- 
dritzki. 
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Tablc 1. All-sky Photometric Magnitudes 



Star*^ 


B V Rr Ir 




XO-1 


11.85 11.19 10.806 10.43 




1 


13.85 12.85 12.282 11.76 




2 


16.31 14.93 14.139 13.41 




3 


12.36 11.54 11.046 10.56 




4 


14.31 13.65 13.313 12.93 




5 


9.99 9.54 9.261 8.97 




6 


12.42 11.29 10.680 10.11 




7 


11.85 10.86 10.330 9.81 




8 


15.08 14.28 13.803 13.28 




a) St 


ars are identified in Figure 1. 




Table 2. The Star XO-1 




Parameter 


Value 


Reference 


RA (J2000.0) 


16'*02'"11^84 


a,b 


Dec (J2000.0) 


+28° 10' 10". 4 


a,b 


V,Vt 


11.19±0.03,11.323 


c,b 


B-V,Bt-Vt 


0.66±0.05,0.692 


c,b 


V - Rc 


0.38±0.04 


c 


Rc-Ic 


0.38±0.04, 


c 


J 


9.939 


A 
Ll 


J-H 


0.338 


J 

a 


H-K 


0.074 


A 
Q 


Spectral Type 


GIV 


C 


d 


200 ± 20 pc 


C 


{l^a,IJ'd) (- 


-19.7 ± 2.0, 15.0 ± 1.9) mas yr'^ 


b 


GSC 


02041-01657 


a 




References: 






a) SIMBAD 




b) Tycho-2 Catalogue, Hog et al (2000) 






c) this work 




d) 2MASS, Skrutskie e al. (2006) 
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Table 3. The Planet XO-lb 



Parameter 


Value 


Notes 


P 


3.941534±0.000027 d 




tc 


2453808.9170±0.0011 (HJD) 




K 


116±9 m s-^ 




Rp/Rg 


(1.30±0.04) X Rj/R© 


a 


a 


0.0488±0.0005 A.U. 


b 


i 


87.7±1.2 deg 


b,c 


Mp 


0.90±0.07 Mj 


b,d 


Rp 


1.30±0.11 Rj 


a,b,c 


Notes: 



a) Rj = 71492 km, i.e. the equatorial radius of Jupiter 

b) for = 1.0±0.03 

c) for R^ = l.OibO.08 R© 
d) Mj = 1.8988e27 kg 



Table 4. Radial Velocity Shifts 



Julian Date Radial Velocity 
Shift [m s-i] 


Uncertainty 
(1 a) [m s-i] 


Telescope 


2453770.947 


-99 


63 


HJS 2.7-m 


2453771.945 


93 


70 


HJS 


2453773.924 


-107 


63 


HJS 


2453775.936 


122 


64 


HJS 


2453779.941 


87 


15 


HET 11-m 


2453796.894 


41 


17 


HET 


2453799.884 


112 


14 


HET 


2453801.887 


-111 


15 


HET 


2453815.848 


128 


21 


HET 


2453824.815 


-34 


18 


HET 
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Table 5. Results of the SME Analysis 



Parameter 


Mean 


Precision 






(1 a) 


TeffiK] 


5750 


13 


\ogg [cm s~^] 


4.53 


0.065 


V sin i [km s"-*^] 


1.11 


0.67 


[M/H] 


0.058 


0.040 


[Na/H] 


-0.051 


0.041 


[Si/H] 


-0.043 


0.038 


[Ti/H] 


-0.000 


0.040 


[Fe/H] 


0.015 


0.040 


[Ni/H] 


-0.060 


0.021 


[Si/Fe] 


-0.058 


0.056 



Table 6. Spectroscopically Derived Stellar parameters 



Parameter 


@ 180 pc 


@ 200 pc 


@ 220 pc 




0.98 


0.97 


0.97 


Mass [Mq] 


1.01 


1.00 


1.00 




1.03 


1.02 


1.02 




0.92 


0.96 


1.04 


Radius [R©] 


0.95 


1.00 


1.09 




0.98 


1.06 


1.15 




4.44 


4.38 


4.31 


Log(g) [cm s-2] 


4.49 


4.43 


4.35 




4.52 


4.48 


4.40 




0.65 


2.48 


5.54 


Age [Gyr] 


2.30 


4.74 


6.93 




4.42 


6.63 


8.01 



For each parameter, the middle row is the maximum hkehhood value, and the 
values in the rows above and below span the 68% likelihood of the probability 
distributions (cf. Figure 7). The three columns correspond to three assumed 

distances for XO-1. 
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Fig. 1. — XO-1 is centered, indicated by the hash mark. Stars from Table 1 are circled. 
North is up; East to the left. The DSS image, digitized from a POSSII-F plate with a IIIaF 
emulsion and an RG610 filter, subtends 19' of declination. 



0.96 I I I 

-2-10 1 2 

Time from mid — transit [days] 

Fig. 2. — More than 3000 individual observations of XO-1 by the two XO cameras over 
two seasons 2004 and 2005 are shown wrapped and phased according to the ephemeris of 
Equation 1 and averaged in 0.01-day bins (hne). Prom these data we identified the star as a 
candidate for more refined photometry with other telescopes at epochs of expected transits 
(Figure 3). 
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-0.05 0.00 0.05 

Time from mid — transit [days] 



Fig. 3. — Time series photometry of XO-1 from 2005, with dates, observers, and fihers 
indicated. The observations have been averaged in 0.006-day bins. The figure is in color in 
the electronic edition. 
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-0.02 -0.01 0.00 0.01 0.02 

Orbital Phase 

Fig. 4. — Time series photometry of XO-1 from the Mar 14, 2006 transit plotted on the 
same scale as Figure 3. We fit the R-band data [diamonds (0.3-m telescope; 0.006 day 
bins) and triangles (1.8-m telescope; 0.002-day bins)] to determine the system parameters 
for XO-1. Although the shape of the R-band light curve provides the better constraint on the 
inclination, the V-band (squares) and B-band (crosses) data (both 1.8-m telescope) appear 
to be slightly deeper and narrower than the R-band, as expected for a near-central {i ~ 90°) 
transit of a limb-darkened star. The solid lines are models appropriate for each filter (color 
coded in the electronic edition as B=blue, V=green, R=red) with = 1.00 M©, = 1.00 
Rq, log g = 4.53, T^fj = 5750 K, = 1.30 Rj, i = 87.7°. 
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200 F 




^ 3770 3780 3790 3800 3810 3820 

Julian Date - 2450000 




0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.5 0.0 0.5 1.0 

Phase sin(Phase) 

Fig. 5. — a) The radial velocity of XO-1 oscillates sinusoidally with a semi-amplitude K = 
116±9 m s~^, implying XO-lb's mass is 0.90±0.07 Mj. b) The period and phase of the radial 
velocities were fixed at values determined by the transits. The mean stellar radial velocity 
with respect to the solar system's barycenter has been subtracted. In order to determine K, 
we used only the HET data cahbrated with an iodine absorption cell (filled circles). The 
data obtained with the 2.7-m HJS telescope and calibrated with ThAr emission fines (open 
squares) are for illustration only. We subtracted mean radial velocities separately from the 
HET and HJS data. Because the data (filled circles and open squares) are from independent 
observers, telescopes, calibration sources, and analyses, we are confident that the radial 
velocity oscillation is genuine, c) In this representation of the data, a circular orbit yields a 
straight line of slope — K. 
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1 — o~T~^ 



—1 I i_ 



6022 6024 
Wavelength (A) 



— 1 1 1 1 — 



— I 1 1 1 1 1 1 1 1 — 

V Fe 



6040 
Wavelength (A) 




6060 
Wavelength (A) 



Fig. 6. — The mean spectrum of XO-1 as observed (black histogram) and modeled with 
SME (curve, colored blue in the electronic edition) in the region of the Mg b triplet. Labels 
note the elements responsible for the indicated spectral lines. Intermittent line segments 
(tan) beneath the horizontal axis indicate wavelength intervals used to constrain the spec- 
troscopic parameters. Very short and intermittent line segments (black) immediately above 
the spectrum indicate wavelength intervals used to constrain the continuum fit. 
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£ 3 
1 2 

CO 

o 

Ql 1 







6.97 1.00 1.02'm. 




= 2 

CO 

£ 1 







0.95 1.00 1.05 

Mass {M,J 

— 1 — 



4.38 4.43 4.48 





0.90 1.00 1.10 1.20 

Radius (R^,,) 

T 



2.48 4.74' 6.63 Gyr 



4.30 4.40 4.50 4.60 

log(g) 




4 8 
Age (Gyr) 



12 



Fig. 7. — Distributions for four stellar parameters derived from the SME analysis (see text) 
for a distance of 200 pc. The values of the mean and limits containing ±34% of the dis- 
tribution from the mean are annotated on the figures and listed in Table 6 along with 
corresponding values for distances of 180 pc and 220 pc. 
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Fig. 8. — The reduced of light curve models rises away from ~1 solar radius, more steeply 
at smaller radii than larger radii because the cosine of the inclination cannot exceed unity 
(cf. Equation 2). Three independent observations of two transits are shown as solid, dotted, 
and dashed lines. The nearly coincident lines correspond to = 0.96, 1.0, and 1.04 Mq 
(colored in the electronic edition). Because the mean stellar density is the critical parameter 
(Equation 2), at a specific x^y the more massive models require larger stellar radii. The 
vertical solid line at r* = IR© is the maximum likelihood value from the SME analysis 
(Table 6); the vertical dotted lines are the corresponding 1-cr limits. 
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Fig. 9. — To this plot of the radii and masses of transiting extrasolar planets in Jovian units 
(cf. footnotes to Table 3), we have added XO-lb (red square) to the latest data for the other 
nine compiled by Charbonneau et al. (2006). Filled symbols are stars brighter than V=12. 
Open circles represent OGLE planets (V~15). Lines of constant mean density are labeled 
in g cm~^. 



